Effect of the introduction of impurities on the 
stability properties of multibreathers at low 
coupling 



J Cuevas, JFR Archilla and FR Romero 

Nonlinear Physics Group of the University of Sevilla 
Departamento de Fi'sica Aplicada I, ETSI Informatica 
Avda Reina Mercedes s/n, 41012 Sevilla, Spain 



Abstract. Using a theorem dubbed the Multibreather Stabiliy Theorem 
[Physica D 180 (2003) 235-255] we have obtained the stability properties of 
multibreathers in systems of coupled oscillators with on-site potentials, with 
an inhomogeneity. Analytical results are obtained for 2-site, 3-site breathers, 
multibreathers, phonobreathers and dark breathers. The inhomogeneity is 
considered both at the on-site potential and at the coupling terms. All the results 
have been checked numerically with excellent agreement. The main conclusion is 
that the introduction of a impurity does not alter the stability properties. 



PACS numbers: 63.20.Pw, 63.20.Ry, 63.50.+X, 66.90.+r. 
E-mail: jcuevas@us.es 



Submitted to: Nonlinearity 
1. Introduction 

Discrete breathers (DBs) are periodic, localized solutions that appear in discrete 
lattices of nonlinear oscillators. For Klein- Gordon lattices, i.e, lattices with an on- 
site potential, the conditions for their existence at low coupling are very weak. They 
are based on properties of the lattice at the anti continuous limit, that is, the same 
lattice with zero coupling and, then, with the oscillators isolated. These conditions 
can be expressed in simple words as: a) For a given frequency Wb the isolated oscillator 
is truly nonlinear (d-E/dwb ^ 0, with E the energy of the isolated oscillator), b) The 
breather does not resonate with the phonons (nwb ^ ui , with n any positive integer 
and lo the frequency of the isolated oscillators . The stability of one-site breathers 
has been proven in |U 03 E] , whereas the stability of two-sites breathers was proven in 
Ref. 0. Recently, some of the authors developed the Multibreathers Stability Theorem 
(MST) , which provides with a method for obtaining the stability properties of any 
multibreather at low-coupling and applied it to homogeneous lattices. Although the 
MST applies to any Klein-Gordon lattice, it depends on the analytical or numerical 
calculation of some magnitudes J nm : This calculation becomes less simple and the 
number of magnitudes larger as the system becomes more complicated, involving, for 
example, different linear frequencies or initial phases for the oscillators. 
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In this paper we address the application of the MST to Klein-Gordon lattices 
with an impurity. This impurity can be modelled at the on-site potential or at the 
coupling. The implementation of it at the masses is equivalent to at the on-site 
potential. The necessary magnitudes are calculated and the MST applied to 2-sites and 
3-site breathers, multibreathers, phonobreathers and dark breathers. The theoretical 
results are also checked numerically to confirm the validity of the calculations with 
excellent results. 

The paper is organized in the following form: first we describe the details of 
the model; second we recall the Multibreathers Stability Theorem, and give some 
details of its application to the case of a lattice with an impurity. Subsequent sections 
particularize it for two and three-sites breathers, phonobreathers, multibreathers and 
dark breathers. Details of analytical calculations are relegated to the first appendix. 
In the second appendix, a conjecture made in Ref. [H] is proven. It is that the only 
multibreathers that can be stable are the ones with all the oscillator in phase or out 
of phase. 

2. The model 

In this paper, we consider Klein-Gordon chains with linear nearest-neighbours 
coupling whose dynamical equations are of the form: 

N 

u n + V'(u n ) + e C nm u n =0 b = 1 JV (1) 

m— 1 

where the variables u n are the displacements with respect to the equilibrium positions, 
V(u n ) is the on-site potential, N is the number of oscillators, C is a coupling matrix 
which includes the boundary conditions, and e is the coupling parameter. 

The impurity can be introduced in the system through an inhomogeneity at the 
on-site potential or the coupling matrix, in a similar fashion as it was done in [7]- 

The linear frequencies of the isolated oscillators are w n = (V£ (0)) 1 / 2 . If the 
inhomogeneity is implemented at the potential of the no-th particle, we denote Vb = V n 
and w — uj n = (^"(O)) 1 / 2 Vn ^= n . Wc can write o; 2 = (1 + ad n>n0 )uj o , where the 
inhomogeneity parameter a takes values in (— l,oo), a — corresponding to the 
homogeneous case. 

For a system with nearest-neighbour interaction, the dynamical equations Eq. |JT)| 
can be written 0: 

il n + V'(u n ) + e[k„-i <n (u n - u n -x) + fc„,„+i (u n - u n+ i)} = 0.(2) 

If the system is homogeneous k n _ rn = 1, \/n, m, except at to = 1 or n — N for a system 
with N oscillators. With periodic boundary conditions, the index n is cyclical so as 
~ N and N + 1 ~ 1 and k 0y i — fcjv,jv+i = 1- With free ends, fco,i = fe/v,Ar+i = 0. 
With fixed ends, two extra oscillators have uo = un+i — and fco,i = fcjv,jv+i = 1- 

Therefore, the coupling matrix, C° for a finite system with N oscillators with 
fixed ends (at the oscillators n = and n = N + 1), has elements: 

{2 if n = to 
-1 if|n-m| = l (3) 
otherwise. 

With free ends, C\\ — Cnn = 1, and with periodic boundary conditions, C\n — 
Cni = — 1. 
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We model the impurity at a site n distinct from the lattice boundaries, by 
changing the constants k n0 n0+ i and fc„o-i,nO to (1 + f3/2), with (3 being a parameter 
which takes its values in (—2, oo), so as the homogeneous case is recovered when j3 = 0. 
f3 < means weaker constants than the homogeneous coupling ones and f3 > the 
opposite. 

Therefore, the coupling matrix C nm becomes: 

([3 + 2) if n = m = n 

. — ((3 + 2)/2 if \n — m\ = 1 and n = no or m = no , 

nm = < ^ (/3 + 4)/2 ifn = m = n ±l () 

C° otherwise, 



3. Stability theorem 

In this section we recall some results of the multibreather stability theorems established 
in Ref. 6 . For more details about this theory, the reader is referred to that reference. 

The stability of a multibreather at low coupling can be established through three 
different parameters: the sign of the coupling constant e, the softness/hardness of the 
on-site potential and the sign of the eigenvalues {A^} of the perturbation matrix Q, 
to be described below. Considering a positive coupling constant, the multibreather is 
stable for a soft potential if all the eigenvalues of Q are non-negative (Q is positive 
semi-dehnite); if the potential is soft, the multibreather is stable if all the eigenvalues 
{Xi} are negative except for a zero one (Q is negative semi-definite). If the coupling 
constant is negative, the results are reversed. Furthermore, if Q is non-definite (i.e. 
there are some positive and some negative eigenvalues), the multibreather is unstable 
independently on the sign of e and the softness/hardness of the potential. 

We can summarize the stability properties in the following way. Let S = 1 mean 
stability and S = — 1 instability, H = 1 correspond to a hard on-site potential, and 
H — — 1 to a soft one, and define the sign(Q) — 1 if all the eigenvalues of Q but a zero 
one are positive and sign(Q) = — 1 if they are negative except for the zero one, other 
cases being indefinite. Then, S — —H x sign(Q) x sign(e). 

It is also important to take into account that there is always a zero eigenvalue 
due to a global phase mode. If there were more than one zero eigenvalue, the stability 
theorem can only predict the instability in the case that there existed eigenvalues of 
different sign, but not the stability, as the 0-eigenvalue is degenerate. 

The perturbation matrix Q is obtained in the following way: suppose that at the 
anticontinuous limit there arep excited oscillators and N — p oscillators at rest. Let us 
construct the modified coupling matrix C by suppressing in C the rows and columns 
corresponding to the N — p oscillators at rest, and redefine C = C. The dimension of 
C is p x p, its nondiagonal elements being j^j: 



T/2 



n m (5) 



t^n J —T/2 

with T — 2t{/l)\>. The diagonal elements are: 

Qnn ^ ^ Qnm • (6) 

V m^n 

where u n are the solutions of an isolated oscillator submitted to the potential V(u n ), 
i.e. the solutions of the equations: 

u n + V'(u n ) - 0, (7) 
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and fi n is denned as: 

/in = v y "" d< (8) 

In this paper we limit ourselves to time-reversible solutions, i.e., the excited 
oscillators at the anticontinuous limit can only have phase or 7r, or, in other words, 
if Unit) is a time-reversible solution of Eq. 10, i.e., u n (t) = u n (—t), then u n (t + T/2) 
is the only possible time-reversible solution apart from u n (t). Therefore, the state of 
the system at the anticontinuous limit can be described by a code a — (ci, . . . ctjv), 
where o~ n = 0, +1, —1 means that the corresponding oscillator is at rest, with phase 0, 
or with phase tt, respectively, at t = [2]. 

Let J nm and J' nm be defined as: 

ii n (t)u m (t)dt, J' nm = [ u n (t)u m (t + T/2)dt, (9) 

-T/2 J-T/2 

and the parameters i] nm , j nm and ip nm as: 



?7nm = — r= > 0) 7nm = ~ — ; > 0, ^nm = ]l ~J > 0, (10) 

V ^ rm mm v ^ nn ^ mm * nn 

Then, the matrix elements Q n m can be written as: 



TJnmCnm if (T n (T m — 1 
^nmCnm if (T n (T m — 



(ii) 



(12) 

V m^n 

where we have taken into account that /i„ = J„ n . We have calculated analytically 
these parameters for several potentials in |Appcndix A| although only the case of the 
Morse on-site potential leads to relatively simple expressions. 

Some properties of the parameters i] nm , 7„ m and tp nm can be easily deduced. In 
the case of symmetric potentials, J nm = —J nm , as u n (t) = —u n (t + T/2), and, in 
consequence, j nm = i] nrn . Other property is that, for the homogeneous case, r/ nm = 1, 
<£mn — 1 and jnm = 70, where 70 is the symmetry parameter in that case (which is 
the unity for a symmetric potential). 

In order to simplify the notation, as we are considering a single impurity, we 
define the parameters 77 = rj n h, 7 = Jnhi V 3 = fnhi where the index n indicates the 
impurity site and h another one. 



4. 2-site breathers with an impurity 

A 2-site breather is a breather obtained from the anticontinuous limit when two 
contiguous sites are excited. The boundary conditions are irrelevant as long as the 
sites are not close to the boundaries. Without loss of generality, let us rename the 
two indexes so as the impurity is located at n = 1, and the corresponding oscillator 
has phase zero (cti = 1). The other oscillator is located at n — 2. When dealing with 
the code a it is enough to consider the codes of these two oscillators, as the other are 
zero, thus, a — (01,02). Two cases are possible, in-phase oscillators, a = (1, 1), and 
out-of-phase oscillators, a = (1, —1). 
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4-1- Inhomogeneity at the on-site potential 



First of all, we consider the in-phase case. Taking into account Eqs. I|1UI12[) and the 
notation introduced in Section [3 the matrix Q is given by: 



Q 



w -v 



with eigenvalues {A , Ai} = {0, r)((p + 1/V)}. Then, Ai > 0. 
For the out-of-phase case, the perturbation matrix is: 



Q = 



-if 

i 



i 

-j/ip 



(13) 



(14) 



with eigenvalues {A ,Ai} = {0,-7(1,9 + 1/V)}- Thus, Ai < 0. 

If the system is homogeneous {Ao, Ai} = {0, 2} and {Ao, Ai} = {0, —270} for the 
in-phase and out-of-phase cases, respectively [6]. In consequence, the sign of Ai does 
not change when an impurity is introduced in the on-site potential, and the stability 
properties are not altered. 

Figs, n and [21 show the dependence of E = e\i with respect to e calculated 
numerically and analytically (see |Appcndix A| ). The agreement between both 
calculations is excellent for fairly large values of the coupling constant. 



4-2. Inhomogeneity at the coupling constant 
For the a = (1, 1) 2-site breather: 



(3 + 2 



1 

-1 



-1 
1 



(15) 



where Qo is the perturbation matrix of the homogeneous chain. Thus, Ai = /3+2. The 
relation Q = (/? + 2)/2Qo also holds for the code a — (1, —1), and, in consequence, 
Xi = —{13 + 2)/27q in this case. As (3 > —2, the conclusion is that the stability 
properties are the same as in the homogeneous system. 

Fig. compares this result with the numerically obtained. 



5. 3-site breathers with an impurity 



They consist of breathers derived from three contiguous, excited oscillators at the 
anticontinuous limit. There are seven non-equivalent possibilities, taking into account 
the position of the impurity and the phases of the oscillators. Let us rename the 
excited sites as 1,2,3, suppose that 02 = 1, and denote with boldface the position of 
the impurity. Then the three-index codes corresponding to the three excited oscillators 
are (1,1,1), (1,1,1), (-1,1,-1), (-1,1,-1), (1,1,-1), (1,1,-1), (1,1,-1)- 



5.1. Inhomogeneity at the on-site potential 

Let us suppose that the a = ((71,1,(73). The in-phase breather corresponds to 
(7i = (73 = 1, the out-of-phase to <7i = 03 = — 1 and another one to (71(73 = — 1. 
The following table shows the eigenvalues Ai and A2 for the different configurations 
(Ao = always), the impurity site is indicated through a bold font, and, in order to 
simplify the expressions, we define = ip + 1/tp: 
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a=-0.2 a=0.2 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.01 0.02 0.03 0.04 0.05 0.06 0.07 




"0 0.01 0.02 0.03 0.04 0.05 0.06 ' 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



e e 

Figure 1. Dependence of E = e\\ with respect to e for two different values of 
a. (a,b) correspond to in-phase 2-site breathers and (c,d) to out-of-phase 2-site 
breathers. All the solutions stand for a Morse potential and u\, = 0.8. Dots 
represent numerical solutions and lines analytical ones. 



Code 


Ai 


A 2 


1 1 1 


v/<p 


v(f + 4>) 


111 


(2 + n<t> + y/4- 4r]<t> + »? 2 2 )/2 


(2 + n <$> + s/4 - A V 4> + V 2 <f> 2 )/2 


-1 1 -1 


-7/V 


-j(>p + 4>) 


-1 1 -1 


-70 - [70 + \f{l4> + 27o) 2 - 477o(¥> + <W]/2 


-70 - hfP - V (70 + 2 7o) 2 - 477o(<^ + </>)]/ 2 


1 1 -1 


{rj<t> + vV</> 2 + ^oviv + 4>))/2 


(ri<f> - vV</> 2 + 4r/ r)((p + <t>))/2 


1 1 -1 


((V ~ 7)* + V(V + 7?<t> 2 ~ 4r m > 2 )/2 


((V ~ 7)0 - v 7 ^ + 7) 2 2 " 4777^ 2 )/2 


1 1 -1 


(-(70 - 2) + v / 7 2 2 + 4 + 47^)/2 


(-(70 - 2) - v / 7 2 </» 2 + 4 + 4 7V )/2 



70 being the symmetry coefficient for the homogeneous case. It can be deduced that 
Ai > and A2 > for the in-phase solution, Ai < and A2 < for the out-of-phase 
breather and Ai > and A2 < otherwise. This result coincides with the homogeneous 
case. 

Figs. 01 and 0] show the dependence of E\ = eAi and E 2 = eA 2 with respect to e 
calculated numerically and analytically. The agreement between both descriptions is 
again very good even for fairly large values of the coupling constant. 
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(c) 
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-0.1 
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Figure 2. Dependence of E = e\\ with respect to e for two different values of 
a. (a,b) correspond to in-phase 2-site breathers and (c,d) to out-of-phase 2-site 
breathers. All the solutions stand for a hard <j> 4 potential and u)^ = 1.3. Points 
represent numerical solutions and lines analytical ones. 



5.2. Inhomogeneity at the coupling constant 

If the impurity is located at no = 2, the relation Q — (/? + 2)/2Qo is fulfilled again. 
The only way to obtain a different relation is to suppose the impurity at an edge. The 
following table summarizes the different cases: 



Code 


Ai 


A 2 


1 1 1 


1 + 0/2 


3(1 + 0/2) 


111 


{13 + 4 + + 4)/2 


(0 + 4 - V0 2 + 4)/2 


-1 1 -1 


-70(1 + ^/2) 


-370(1 +0/2) 


-1 1 -1 


-70 (0 + 4+ ^0 2 +4)/2 


-7„(/3 + 4+ V0 2 +4)/2 


1 1 -1 


1 - 70 + (0 + sj 7,)0 2 + 2 7 () (270 + 1)0 + 4( 7 2 + 70 + l))/2 


1 - 70 + (0 - y/ 7,)0 2 + 270(270 + 1)0 + 4( 7 2 + 70 + l))/2 


1 1 -1 


(1 + 0/2)(1- 7o + /yo +70 + 1) 


(l + 0/2)(l- 7 o + v /7o +70 + 1) 


1 1 -1 


1 - 70 + (0 + ^700" + 270(270 + 1)0 + 4( 7 2 + 70 + l))/2 


1 - 70 + (0 - 7 7o0 2 + 270(270 + 1)0 + 4( 7 2 + 70 + l))/2 
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The signs of the eigenvalues are the same as in the on-site inhomogeneity case. Fig. [S] 
compares these results with the numerically obtained ones. 

6. Multibreathers and phonobreathers with an impurity 

A multibreather is obtained from the anticontinuous limit when a number of 
contiguous oscillators are excited. If all the oscillators are excited it is called a 
phonobreather. Provided that a multibreather does not include the oscillators at 
the borders, the modified coupling matrix of a N-site multibreather is identical to 
the one of a phonobreather with fixed ends. Therefore, we can study both cases 
simultaneously. We cannot obtain analytical expressions for the eigenvalues of the 
perturbation matrixes, but we have been able to demonstrate their stability properties. 

When the number of oscillator of a multibreather increases, the degree of the 
characteristic polynomial also increases and an analytical evaluation of its roots is 
not possible. However, in the case of a homogenous system, the eigenvalues can 
be calculated because the eigenvectors equation of the system is equivalent to the 
dynamical equations of a linear lattice of oscillators, as shown in Although 
this analogy might be suggested in our case, the only information available for the 
eigenvalues would be the dependence on the parameters of a localized eigenvalues and 
only in an infinite lattice. 

Thus, we can only obtain qualitative results about the spectrum of the matrix Q 
for a phonobreather (or N-site multibreather). In particular, we show that the stability 
properties of the system with an impurity are the same as an homogeneous system. 
In order to do that, we will make use of some congruence properties of symmetric 
matrices [S]- 

Two symmetric matrices A and B are congruent if they can be transformed each 
other through elementary transformations: B = PAP 1 . The inertia of a matrix 
is defined as In(A) = {i+(A),i-(A),io(A)} where i+, i-, io denote, respectively, 
the number of positive, negative and zero eigenvalues of A. Sylvester's inertia law 
establishes that the inertia of two congruent matrices are the same [H] . In consequence, 
it is enough to diagonalizc a matrix using elementary transformations to obtain its 
inertia. The diagonal matrix has the structure called first canonical form: 



We analyze below the stability properties of phonobreathers with different 
boundary conditions and an inhomogeneity at the on-site potential. Let us recall 
that, in the homogeneous case, the inertia of Q is, for an in-phase breather, In(Q) = 
{N — 1,0,1} and, for an out-of-phase (staggered) breather, In(Q) = {0,iV — 1,1}, 
with N being the number of system particles. We will show that the inertia does not 
change when an impurity is introduced. 

Furthermore, the only multibreathers that can be stable are those that vibrate 
in-phase or staggered, as demonstrated in |Appcndix B| 

6.1. Multibreathers and phonobreathers with free/fixed ends boundary conditions 

If the impurity is supposed to be located at n = c, the elements of the perturbation 
matrix Q F for an in-phase multibreather, can be written as: 




(16) 
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(a) (b) 

a=-0.1 a=0.1 




£ £ 

Figure 3. Dependence of E = eXi and E = e\2 with respect to e for two different 
values of a. (a,b) correspond to in-phase 3-site breathers, (c,d) to out-of-phasc 
3-site breathers with the impurity at the center and (e,f) to the same as before 
but with the impurity at the edge. All the solutions stand for a Morse potential 
and ojt, = 0.8. Points represent numerical solutions and lines analytical ones. 
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Figure 4. Dependence of E = e\\ and E = e\2 with respect to e for two different 
values of a. (a,b) correspond to in-phase 3-site breathers, (c,d) to out-of-phasc 
3-site breathers with the impurity at the center and (e,f) to the same as before 
but with the impurity at the edge. All the solutions stand for a hard <fr 4 potential 
and cjjj = 1.3. Points represent numerical solutions and lines analytical ones. 
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Figure 5. Dependence of E = eAi for a 2-site and a 3-site breather with an 
impurity at the coupling located at the edge. All the solutions stand for a Morse 
potential and ui^ = 0.8. Points correspond to numerical solutions and lines to 
analytical ones. 



Qnm = 2<5 n ,m[l + (w - 1)<W + (v/<P ~ l)/2<5|„_ c |,i - (£„,! + 5 n<N )/2] - 

-<5n,m+l[l — (V + l)(^n,c + <W+l)] - £n,m-l[l - (j? + l)(6 n<c + 6 n ,c-l)]- (17) 

This matrix is transformed into a diagonal matrix through D = PQP 1 , where 

Pnm = [8[2i],0 + Sn,m][l + ^],0^n,c-l(V <p/V ~ 1) + ^m,c<%],0 / V + 

+ <y„, c -i((vW - i)*[?],o + (i/v^j - i)<W)L (is) 

where [•] denotes the integer part. D nm = S„ tm (l — 5 n ,iv)- I n consequence, 
In^^ ) = {TV — 1, 0, 1} and all of the eigenvalues are positive except for one which is 
zero. 

The perturbation matrix for the out-of-phase breather is given by: 

Qnm = -270<W1 + (7V - l)*»,c + (7/V ~ l)/2*|n-o|,l ~ («M + <W)/2] - 
+70<^n,m+l[l - (70 + l)(5„, c + 5 n ,c+l)] + 70<Wn-l[l ~ (7/70 + 1)(<W + ^n.c-l)], (19) 

which transform into the diagonal matrix D nm = — S n<m (l — 5 n ,jv) through the 
transformation matrix 

-Pnm = [8[M],0 + ^n,m][l + %],0^,e-l(\/WT _ 1) + &m,cS[a],0 / <P + 

+ S n ,c-i((VWl - l)%],o + (VV^7 - 1)<W)]. (20) 

The inertia of the stability matrix is given by ln(Q Fn ) = {0, N — 1, 1} and, in 
consequence, it has all its eigenvalues negatives except for the null one. 

Figure [B] shows the numerical values of the eigenvalues and confirms the results 
stated previously. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

a 



Figure 6. Dependence of the eigenvalues of the perturbation matrix with respect 
to the impurity parameter a for a phonobreather with free/fixed ends boundary 
condition, (a) and (b) corresponds to the Morse potential and = 0.8 while (c) 
to a hard </> 4 potential, (a) and (c) corresponds to a in-phasc breather and (b) to 
a staggered one. 

6.2. Periodic boundary conditions 

The diagonalization of Q is now very difficult. We will make use instead of a theorem 
due to Weyl [B] to obtain the inertia of the matrix. This theorem establishes that, if 
the eigenvalues of a matrix are arranged in the following order: Ai < A2 < . . . < A^r 
then, the k-th eigenvalue of the sum of two matrices holds: 

X k (A) + Ai(B) < X k (A + B)< X k (A) + X N (B). (21) 

This theorem can be particularized for our case if we denote as Q P0 and Q Pt! to 
the periodic boundary conditions matrices and define Q p = Q F + A, where Q F is the 
perturbation matrix with free boundary conditions. By making A = Q F and B = A 
we can apply Weyl's theorem. 

For a phonobreather in-phase, 

= 2(5„, m [l + (rjip - l)5 n . c + (r)/cp - l)/2(S|„_ c | jl ] - {5 n ^8 m . N + S ntN S m< i) - 

S n<m+1 [l - (?7 + l)(£„, c + <W+l)] ~ <Wn-l[l - ( r l + 1)(<W + Sn,c-l)], (22) 
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and, in consequence, 

A„ >m = 5 n ,m\5 n ,l + 5 n ,N) — (o~n,lO~m,N + <5n,AT(5 m: i). (23) 

The spectrum of A is given by spec(A) = {0, 2} where Ai(A) = has multiplicity 
N — 1 and Ajv(A) = 2 has multiplicity 1. Thus, particularizing 1|21JI for k = 1: 

Ai(Q F0 ) + Ai(A) < Ai(Q P0 ) < A!(Q F0 ) + X N (A). (24) 

As Ai(Q F0 ) = 0, Ai(Q F0 ) € [0,2]. If k = 2 is taken in (ZD, A 2 (Q F0 ) > 
A 2 (Q F0 ) + Ai(A). As A 2 (Q F0 ) > then A 2 (Q F0 ) > 0. In consequence, Q P0 has 
N — 1 positive eigenvalues and one which is greater or equal to zero. 

The steps for an out-of-phase multibreather are similar to the in-phase case, 
except that k = N is taken in 1|21|1 . Thus, it is straightforward to show that Q P7V has 
all its eigenvalues negatives except for one of them which is smaller or equal to zero. 

Figure0shows the numerical values of the eigenvalues and confirms the previously 
stated results. 

Let us remark that for N odd there appear, in a similar fashion to the 
homogeneous case, parity instabilities due to the breaking of the breather pattern. 
Figure |H1 shows the dependence of the eigenvalues of Q with respect to a in this case. 
It can be observed the existence of a positive eigenvalue of constant value. It was 
demonstrated in ^U] that in a homogeneous lattice, this eigenvalue is ij /3 in an 
infinite lattice. This value does not change when an impurity is introduced as it is 
due to a border effect. 

6.3. Inhomogeneity at the coupling constant 

In this case, the perturbation matrix is equivalent to the one corresponding to an 
inhomogeneity implemented at the on-site potential with the changes <p = 1 and 
?7 = 7 = 1 + (3/2. In consequence, the results are qualitatively the same in both cases. 

7. Dark breathers with an impurity 

A particular case of multibreathers in which all the particles of the lattice are excited 
except one is called a dark breather 1 1 II I12| . In the case that there are p adjacent 
sites with no vibration, we deal with a p-site dark breather. For the sake of simplicity 
we will consider 1-site dark breathers as the generalization to p-site dark breathers is 
straightforward. 

It is worth noticing that if the impurity is located at the dark site, the perturbation 
matrix turns into the homogeneous one. In consequence, the only way to introduce 
a modification in the perturbation matrix is to suppose the impurity adjacent to the 
dark site. 

For an in-phase breather, the perturbation matrix is given by (notice that the 
c + 1 row and column are removed) : 

Qnm = 2<W1 + (W - l)*n,c + (v/<P - l)/2*n,c-l + Q> - l)(S„,l + 5 n , N )/2] - 

-S n , m+ i[l - (r) + 1)£„, c _i] - <5„, m -i[l - (?7 + 1)<5„ jC _i] - b(Si tN + Si tN ), (25) 

where b — for fixed/free ends and 6 = 1 for periodic boundary conditions. If 
6 = the perturbation matrix is block diagonal, being the lower block the matrix of 
a phonobreather with fixed/free ends boundary conditions and N — (c + 1) particles. 
The upper block is a matrix Q\ with positive eigenvalues except for a null one. In 





Figure 7. Dependence of the eigenvalues of the perturbation matrix with 
respect to the impurity parameter a for a phonobreather with periodic boundary 
condition, (a) and (b) corresponds to Morse potential and ui^ = 0.8 while (c) to 
a hard (j> 4 potential, (a) and (c) corresponds to a in-phase breather and (b) to a 
staggered one. 



consequence, the inertia of the perturbation matrix for b = is {N — 2,0,2}. The 
existence of two zero eigenvalues makes it impossible to determine the stability through 
the theorem as, to first order of the perturbation, the degeneration is not raised. This 
result was also obtained in the homogeneous case. 

For 6=1, i.e. periodic boundary conditions, the perturbation matrix can be 
transformed into Qi through row/columns interchange. This transformation lets the 
characteristic polynomial invariant. In consequence, the matrix has N — 1 positive 
eigenvalues and a null one. 

The perturbation matrix of a staggered dark breather is 

Qnrn = -2-yo6 n , m [l + {~f<p - l)tf n ,c + {l / f ~ l)/2<J n ,c-l + (6 - l)(<fn,l + 6 n , N )/2] - 

- 70<5n,m+l[l - (7/70 + l)5n,c-l] 

- 7o5n,m-l[l - (7/.90 + l)5n,c-l] - &7o(5l,JV + Si,n)- (26) 

It can be shown straightforwardly that In(Q D7T (b = 0)) = {0,7V - 2,2} and 
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0.5 1 1.5 2 0.5 1 1.5 2 

a a 



Figure 8. Dependence of the eigenvalues of the perturbation matrix with respect 
to the impurity parameter a for a Morse potential and uj h = 0.8 in the case 
of staggered (a) phonobreathers and (b) dark breathers with periodic boundary 
conditions. The existence of an isolated positive eigenvalue A = 47o/3 ~ 1.07 is 
observed. 



ln(Q D7T (b = 1)) = {0, N - 1, 1}. In the case of N odd and 6=1, parity instabilities 
also occur (see Figure |SJ. 

Figure |5] shows the eigenvalues for a dark breather with periodic boundary 
conditions. 

8. Conclusions 

The Multibreathers Stability Theorem [S] provides with a method for obtaining the 
perturbation matrix Q from the anticontinuous limit of any multibreather at low 
coupling in Klein-Gordon systems. The stability depends on the signs of the 
eigenvalues of Q, the hardness/softnes of the on-site potential and the sign of the 
coupling parameter (i.e., the attractiveness or repulsiveness of the coupling). In this 
paper we apply the MST to systems with an impurity either at the on-site potential 
or the coupling constants, the first case being equivalent to an impurity at the masses. 

Using analytical methods we have been able to obtain the eigenvalues of Q for 2- 
site and 3-site breathers and compare them with numerical results showing an excellent 
agreement. For larger multibreathers, phonobreathers and dark breathers we have 
obtained analytically the signs of the Q-eigenvalues and checked them numerically. In 
all cases, although the values of the eigenvalues change, the stability properties of the 
system with an impurity coincides with the ones of the homogeneous system. 

The necessity of either an in-phase pattern or an out-of-phase one for the stability 
of a multibreather in the homogenous system, as it was conjectured in Ref. has 
also been demonstrated. 
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O 1 ~ ^ iii 
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

a 



Figure 9. Dependence of the eigenvalues of the perturbation matrix with respect 
to the impurity parameter a for a dark breather with free/fixed ends boundary 
condition, (a) and (b) corresponds to Morse potential and tjy, = 0.8 while (c) to 
a hard <fi 4 potential, (a) and (c) corresponds to a in-phase breather and (b) to a 
staggered one. 



Appendix A. Analytical calculation of parameters 

The aim of this appendix is to calculate analytically the parameters rj nm , j nm , <Pnm 
defined in Section as: 

Jnm „ Jnm - n / J™,m „ f k -t\ 
r) nm = > 0, 7„ m = > 0, (p nm = \ — > 0, (A.l) 

V "iwi^mm V "Tin tJ 'mm V Jnn 

with 

Jnm = [ u n (t)u m (t)dt, J' nm = [ u n (t)u m (t + T/2)dt. (A.2) 

J-T/2 J-T/2 

Analytical calculation of these integrals is rather difficult. In order to overcome 
this difficulty solutions u n (t) are expressed as a Fourier series expansion: 

oo 

U n (t) = Zfj.n + 2 X! Zk - n C0S ( kuJ ^ t), (A. 3) 

fe=l 
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and the integrals turn into series: 

oc oo 

Jnm = 47rw b ^ k 2 z kiTl z kim , J' nm = 47rw b ^(-l) fe fc 2 z fci „z fcjm . (A. 4) 

k=l k=l 

In consequence, the theorem parameters can be expressed as a function of the 
Fourier coefficients of the solution. In the case of the Morse potential, the sum can be 
performed and a closed expression is obtained. 



Appendix A.l. Morse potential 

The Morse potential has the expression V{u n ) = (1 + a n )(exp(— u n ) — l) 2 /2, and the 
orbits of the equation il n + V'(u n ) = are given by 



u n (t) = loj 



1 T \/l - (LOb/Un) 2 COSUJbt 



(A.5) 



The minus sign corresponds to u n (0) > and the plus sign to u n (0) < with 
= V 1 + OL n ■ The Fourier coefficients are: 



Zo, n = log 



2^ 



Zk.r 



Pk I W n - LU h 



k \ UJ n + Ld h 



fc/2 



with pfc = (-l) fe+1 if u„(0) > and p k = -1 if tt n (0) < 0. Thus, 



fc=l 
oo 



J' = 



47TW b ^ 



fc=l 



(0J n + W b )(w m + W b ) 
(C^n - a>b)(^m - ^b) 

(w n + u b )(w m +w b ) 



fe/2 



fc/2 



(A.6) 



(A.7) 



(A. 



The sums are geometric series whose ratios have absolute values smaller than one. In 
consequence, they can be easily summed and the resulting values of J and J' are: 

j» m = 4^ b ; VK-^K,^) (Ag) 

y/ (w n + Wb)(w m + Wb) - V (w tt - Wb)(w m - Wb) 
■^nm = 47TWb- 



^/(w n + w b )(w m + CJh) + y/(uj D - w b )(w m - w b ) 
Then, applying equation l|A.l|) . the theorem parameters are: 

2oj h 



Tfrvm 



Trim 



^/ (w„ + Wb)(u m + CJ b ) - -l^b)(Wm - Wb)' 

2^ 

a/ (W n + W b )K + LOb) + VK -Wb)(w m - U>b) 



U m — Wb 
Wn - ^b 



(A.10) 

(A.11) 
(A.12) 

(A.13) 



Appendix A. 2. </> 4 potential 

This potential is given by l^(u„) = (1 + a Il )ii 2 [ /2 — su^/4, with s = ±1. The sign of 
s allow us to define two different cases: 
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Appendix A. 2.1. Hard <f) A potential, s = +1. The orbit of a hard oscillator is given 
by: 



i n (t) = ±W n i 



2n 2 



LU n t 



2kI (2K(K n ) 



11-2^ V 7T 



(A.14) 



where cn is a Jacobi elliptic function of modulus n n and K(n n ) is the complete elliptic 
integral of the first kind defined as: 

-/ 2 dx 



K(k) = 



\/ \ — n 2 sin 2 x 



The breather frequency Wb is related to the modulus K n through: 

7TW n 



(A.15) 



(A.16) 



2y/l - 2K 2 n K( Kn ) 

The elliptic function can be expanded into a Fourier series and it is obtained 

/ 9 ^ +1 / 2 

v = 0,1,2, (A.17) 



g„ is the Nome and is defined as 



g„ = cxp(- TrK(y/l - «£)/#(«„)). 
Thus, J nm is given by (note that J' = — J as the potential is even): 
87T 3 cj b / q n q m ^2, (2u + \) 2 {q n q m ) u 



Jnm ~ K( Kn )K( Km )\ , 
and, in consequence 



q n q m \ - \£v + i, 

l-2 K 2)(l-2^)^(l + g r +1 



(A.18) 
(A.19) 



(2^+l) 2 (g n g m )^ 



(A.20) 



f,/ (2,+ l)^ V " /(2 



(2^ + 1)C 



9 



2i;+l 



^7, 



Lo m K(n n ) 

LU n K(K m ) 



i/=0 ^ 



(2* + 1)C 



2i/+l 



2 • 



(A.21) 
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Appendix A. 2. 2. Soft </> 4 potential, s = — 1. Now, the orbit of an oscillator is given 
by: 



u n (t) = ±w n . 



2K icd 



UJ„t 



l + Kl 



2«2 f2KU n ) . t 

^rCd I — W b Kn , (A.22) 



1 + «g 



where cd is another elliptic function. Now, the breather frequency co>b and the modulus 
Kjji are related through: 

TTLU 

UJ h = 



Following the Fourier coefficients can be obtained: 



±- 



v+l/2 



2v+l ' 



and, in consequence 

87T 3 Wb 



J, 



q n q m \ - (2^ + l) 2 (g n q m )" 

V (1 + «£)(! + K 2J 2-, (1 „ ^+i )fl „2„+l 



i/=0 ^ 



)^ (i-<zr +i )(i-^ +1 )- 

{2v+\) 2 {q n q m y 



L )(l-< +i ) 



^/(2^ + 2M 

\ Z^l i_ ^+1 

W m K(K n ) 
LO n K(K m ) 



9 



2i;+l 



(i + <4)£(« 

„ \ 1 — Qm 



u=0 



\ !/=0 V 



(2f + 1)<£ 



3i 



2i/+l 



2 ' 



(A.23) 



(A.24) 



(A.25) 



(A.26) 



(A.27) 



Appendix B. Demonstration of the necessity of an in-phase or 
out-of-phase pattern for the stability of multibreathers 

It was shown in Ref. |o] that a multibreather with nearest neighbour coupling and 
whose vibration pattern is neither in-phase nor out-of-phase is unstable independently 
of the hardness/sotftness of the on-site potential and the sign of the coupling constant. 
However, this fact was not demonstrated. 

In this appendix, we will demonstrate the last assertion by making use of a 
consequence of Sylvester's inertia law, called Sylvester's Theorem It establishes 
that, for a positive definite matrix, the principal minors are positive, for a negative 
definite one, the principal minors alternate their signs when the dimension increases by 
one and, finally, a matrix is not definite when none of the last conditions are fulfilled. 

Let us suppose an homogeneous lattice. If the pattern of vector a is different from 
the in-phase or staggered ones, there is always a sequence { — 1, 1,1} or {1,-1,-1} on 
it. It is easy to show that in both cases, the perturbation matrix has the form: 



(B.l) 
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with 

r -70 70 

Q°= 7 l-7o , (B.2) 
1 • 

where • is a number that depends on the phase of the particle adjacent to the 
pattern breaking sequence. As the spectrum of Q is invariant under rows/columns 
exchanges, the block Q° can be placed at Q1.1, and Sylvester's theorem can be applied 
to Q°. The 1st order minor is = Q\ 1 = —70 < 0, and the 2nd order one, 

M^ 2 ) = Qi 1Q22 ~ Q\ 2Q2 1 = ~7o < 0. In consequence, the matrix is not definite 
and there are positive and negative eigenvalues. 

Let us suppose now that an impurity is introduced in the pattern breaking sites. 
The matrices for impurities in the first, second and third sites are (note that we 
only consider the first (2 x 2) submatrix as the result is independent on the third 
row/column), respectively: 
7 

7 -7/V+l 



Qi 



Q°2 



7o 



(B.3) 



-7/V 7 Q° =\ ~ 7 ° 

7 {V- 7)/V J 3 L To v/f - 7o 

and the corresponding minors are: {M^ = —jip,M[ 2 ^ = — "ftp}, {M^ = 
— j/f, = — 777}, {M^ = —jo,M^ = —j r]/(p}. In consequence, the matrices 
are not definite. 
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